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Preface 


Often  the  research  or  military  communities  require  a  quick  estimation  of  neu¬ 
tron  radiation  on  human  body  organs.  A  method  of  estimation  can  help  solve 
problems  ranging  from  amount  of  radiation  shielding  required  in  a  laboratory  to 
distance  required  to  protect  air  crews  from  their  own  nuclear  bomb  drops. 

I  would  like  to  thank  Major  Beller  for  presenting  a  thesis  topic  I  actually 
liked  —  then  letting  me  work  on  it.  I  also  thank  him  for  his  help  and  guidance 
using  the  sometimes  "trying"  MORSE  code.  Special  thanks  to  my  cronies  for  put¬ 
ting  up  with  me  in  times  of  stress  and  taking  me  out  hunting  and  fishing  for  the 
much  needed  breaks.  An  added  thanks  to  the  mutts,  J.D.  and  Magnum,  for 
warming  my  feet  while  I  typed  and  distracting  me  when  I  needed  it.  And  above 
all,  I  thank  my  parents  who  made  this  all  possible,  and  Vicki  who  made  it  all 
bearable. 
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Abstract 

4 

This  report  describes  a  computer  method  of  determining  absorbed  neutron 
dose  to  a  human  phantom.  Modifications  to  the  Oak  Ridge  National  Laboratory 
MORSE  Monte  Carlo  code  result  in  a  code  capable  of  estimating  absorbed  dose  on 
a  human  phantom  in  the  standing  position.  The  phantom  organs  analyzed  are  the 
skin,  bone,  brain,  gastro-intestinal  tract,  and  all  remaining  tissue.  The  organ 
choices  are  limited  to  organs  capable  of  incapacitating  a  human.  The  code  allows 
for  five  different  source  direction  configurations  that  simulate  neutrons,  of  any 
specified  energy  distribution,  incident  on  the  phantom. 

MORSE  analysis  of  a  fission  neutron  spectrum  on  the  phantom  produces 
absorbed  dose  estimates  comparable  with  Japanese  atomic  bomb  survivor  dose 
estimates  by  Scientific  Applications  International  Corporation.  The  analysis  of 
24,001)  source  neutrons  requires  less  than  15  central  processing  unit  minutes  on  a 
VAX  11/780  computer  (VMS  operating  system).  Although  the  code  is  currently 
usable,  additional  phantom  model  orientations,  energy-dependent  quality  factors, 
and  implementation  of  secondary  gamma-ray  dose  estimation  could  greatly 
improve  the  flexibility  and  usefulness.  J  lei.  1  ■  ■  r  ■ 


A  METHOD  FOR  DETERMINING  NEUTRON  DOSE 
TO  A  HUMAN  PHANTOM 


I.  Introduction 


Background 

The  nuclear  age  brought  with  it  many  new  problems  for  both  research  scien¬ 
tists  and  military  strategists  alike.  One  such  problem,  and  the  focus  of  this 
research  effort,  is  determining  the  absorbed  radiation  dose  to  a  human  fiom  neu¬ 
trons  produced  in  reactors,  weapons,  or  simple  neutron  sources.  Effectively  solv¬ 
ing  this  problem  for  organs  capable  of  incapacitation  by  neutrons  could  be  of 
great  interest,  especially  to  military  strategists,  for  the  following  and  many  more 
reasons: 

1.  A  "safe"  distance  between  nuclear  targeted  locations  and  ground  troops 
must  be  known  to  ensure  troop  survivability  in  a  theater  of  conflict. 

2.  Air  crews  delivering  nuclear  weapons  to  many  targets  must  be  able  to 
reach  minimum  "safe"  distances  from  detonation  locations  to  survive 
effects  from  exposure  to  each  successive  detonation. 

3.  The  amount  of  shielding  required  in  above  or  below  ground  shelters 
must  be  known  to  protect  individuals  from  other  than  "safe"  radiation 
environments  so  that  operations  are  sustained  throughout  combat  condi¬ 
tions. 


In  all  three  cases,  "safe"  refers  to  the  maximum  amount  of  radiation  the  average 
human  can  withstand  without  suffering  incapacitation  during  the  length  of  a 
required  duty  (i.e.  a  B-52  bomber  flight). 

The  scientific  community  possesses  many  methods  of  estimating  neutron 
doses  for  specific  cases  such  as  the  Hiroshima  and  Nagasaki  atomic  bomb  victims 


(6:610-*’  19).  Other  studies  (11)  provide  dosimetry  calculations  from  radiation 
effects  correlations.  The  list  of  studies  and  estimations  is  endless.  Unfortunately, 
all  of  the  studies  are  specific  and  time  consuming. 

Problem  and  Scope 

The  problem  is  to  determine  a  fast,  convenient,  and  generic  method  of 
estimating  the  radiation  dose  to  specific  human  organs  from  neutrons  produced  in 
a  nuclear  environment.  The  focus  of  this  thesis  is  twofold.  First,  the  study 
creates  a  human  phantom  capable  of  accurately  representing  the  human  body 
with  organs  of  interest  to  include  skin,  bone,  muscle,  brain,  and  gastro-intestinal 
tract.  Second,  the  study  produces  a  user-ready  version  of  the  Oak  Ridge  National 
Laboratory  (ORNL)  MORSE  computer  code.  MORSE  is  a  Monte-Carlo  method, 
neutron/gamma-ray  transport  code  distributed  by  the  Radiation  Shielding  Infor¬ 
mation  Center  (RSIC)  at  ORNL.  The  final  modified  version  of  MORSE  should  be 
capable  of  determining  a  normalized  absorbed  organ  dose  in  a  human  phantom 
from  any  of  five  source  configurations: 

1.  isotropic  from  a  sphere  surrounding  the  phantom, 

2.  mono-directional  from  the  right  side  of  the  phantom  (+x  direction), 

3.  mono-directional  from  above  the  phantom  (-z  direction), 

4.  mono-directional  from  the  front  of  the  phantom  (+y  direction), 

5.  mono-directional  from  behind  the  phantom  (-y  direction). 

Further,  the  final  version  should  allow  for  a  variety  of  source  energy  spectra  to 
allow  for  various  problems  ranging  from  thermonuclear  or  fission  weapons  to  sim¬ 
ple  laboratory  neutron  sources. 

The  code  modification  is  limited  to  a  standing  phantom  with  five  defined 
organs:  skin.  bone,  brain,  gastro-intestinal  tract,  and  muscle  (includes  the 
remainder  of  the  body  tissue).  Only  these  five  organs  are  modeled  since  the 


primary  interest  is  an  immediately  incapacitating  dose. 


Approach 

The  ultimate  goal  of  this  project  is  reached  by  following  six  major  steps: 

1.  development  of  a  three-dimensional  phantom  detailed  enough  to  provide 
accurate  representation  of  the  human  body  and  the  organs  of  interest. 

2.  preparation  of  the  multi-group  neutron  cross  sections  required  to  model 
the  phantom  organs  and  surrounding  air  volume, 

3.  familiarization  with  the  standard  version  of  the  MORSE  computer  code, 

4.  modification  of  MORSE  subroutines  to  produce  a  phantom-problem 
specific  code, 

5.  analysis  and  verification  of  the  phantom  problem  for  a  fission  weapon 
neutron  spectrum, 

6.  final  preparation  of  the  user-ready  version  of  MORSE. 

The  first  and  second  steps  are  self  explanatory.  The  third  step  is  essentially 
debugging  MORSE  and  running  a  sample  problem  included  with  the  RSIC  code 
package.  The  fourth  step  covers  five  stages  in  the  evolution  of  the  phantom  ver¬ 
sion  of  MORSE  using  various  geometric  shapes  to  test  the  subroutine 
modifications.  The  modifications  are  detailed  in  the  next  chapter. 

The  fifth  step  is  an  exercise  in  producing  results,  testing  the  sensitivity  of  the 
results  to  variable  parameters,  and  validating  the  results. 

The  final  step  is  the  preparation  of  the  user-ready  version  of  MORSE.  The 
user-ready  version  is  the  goal  of  the  thesis  since  it  provides  a  means  of  easily 
determining  normalized  absorbed  organ  doses  for  a  wide  variety  of  problems. 


1-3 


Sequence  of  Presentation 

The  remaining  pages  of  this  thesis  describe  the  stages  of  work  in  detail. 
Chapter  II  describes  the  phantom  and  cross  section  production  followed  by  the 
evolution  of  the  MORSE  subroutines.  Additionally,  Chapter  II  includes  some 
background  about  MORSE.  Chapter  III  contains  the  results  of  this  research  along 
with  a  sample  case  example.  The  final  chapter  (IV)  contains  the  conclusions  of 
this  thesis  and  recommendations  for  further  research. 


II.  Preparation  and  Evolution  of  MORSE  Files 

* 

Overview 

The  MORSE  computer  code  package  is  a  complex,  and  sometimes  confusing, 
menagerie  of  FORTRAN  subroutines  and  input  data  files.  Analyzing  a  given 
problem  may  require  modifications  to  many  of  the  subroutines  and  data  files  --  a 
feat  not  easily  performed  without  understanding  MORSE  and  the  Monte  Carlo 
method. 

This  chapter  details  the  steps  performed  to  modify  MORSE  to  meet  the  goals 
of  this  thesis.  First,  some  background  about  MORSE  is  presented  followed  by 
sections  on  the  preparation  of  the  human  phantom  and  neutron  cross  sections. 

The  last  portion  of  the  chapter  includes  the  10  stages  involved  in  the  evolution  of 
the  final  MORSE  version. 

MORSE  Background 

The  MORSE  (Multi-group  Oak  Ridge  Stochastic  Experiment)  code  is  a 
multi-group  neutron  and  gamma-ray  transport  code  developed  at  ORNL  (13:4.2- 
1).  The  code  is  based  on  Monte  Carlo  methods  for  determining  probable  particle 
interactions  for  a  wide  variety  of  problems  in  three-dimensional  geometry.  Solving 
problems  in  three-dimensional  geometry  is  made  possible  by  the  combinatorial 
geometry  (CG)  module  in  MORSE.  The  CG  module  determines  three-dimensional 
volumes  and  boundaries  from  combinations  of  various  geometric  shapes  as  dis¬ 
cussed  in  the  MORSE  applications  guide  (5:21-28). 

MORSE  uses  Monte  Carlo  probability  distribution  functions  to  determine 
fluence  estimates  from  a  user  specified  sample  population.  MORSE  tracks  the 
entire  population  from  particle  "birth"  (source  production)  through  particle 
"death"  (escape  from  geometry,  variance  reduction  kill,  or  time  kill,  all  discussed 


below)  assigning  a  statistical  weight  at  each  stage  of  the  particle’s  life. 

A  particle  birth  includes  a  starting  location,  energy  group,  direction,  and  sta¬ 
tistical  pathlength.  Along  the  pathlength  a  particle  collides  with  host  nuclei  or 
crosses  boundaries  between  two  geometry  zones.  Following  each  collision,  MORSE 
multiplies  the  particle’s  weight  by  Ea/Et  (£s  and  Et  are  the  scattering  and  total 
macroscopic  cross  section,  respectively)  to  account  for  the  probable  fraction  lost  to 
absorption.  Additionally,  MORSE  determines  a  new  energy  and  direction  follow¬ 
ing  each  collision. 

Particles  crossing  a  boundary  into  a  defined  region  either  undergo  collisions 
or  cross  another  boundary.  Particles  leaving  (escaping  from)  the  defined  boun¬ 
dary  are  "killed"  by  reducing  their  weight  to  zero.  Figure  1  is  a  simplified  flow¬ 
chart  of  a  particle  history  (referred  to  as  "random  walk”). 

Unfortunately,  if  a  particle  does  not  leave  the  geometry  it  continues  to 
"bounce"  around  to  very  low  energies  with  little  or  no  contribution  to  an  overall 
fluence  estimate.  The  many  low  energy  collisions  can  add  substantial  computa¬ 
tional  time  to  the  overall  problem.  On  the  other  hand,  if  particles  escape  from 
the  system  at  high  energies,  smaller  geometry  zones  may  not  receive  enough  track- 
lengths  to  produce  good  statistics.  Time  and  variance  reduction  techniques  are 
available  with  MORSE  to  effectively  solve  both  problems.  Time  kill,  Russian 
roulette,  and  splitting  are  examples  of  such  techniques  and  are  described  below. 

The  time  kill  technique  is  simply  a  particle  age  limit  determined  by  the  user. 
Throughout  the  random  walk  MORSE  keeps  track  of  each  particle’s  age.  When 
the  age  of  a  particle  is  longer  than  the  user-specified  time  kill  value,  the  particle’s 
weight  is  set  to  zero. 

Russian  roulette  is  also  useful  for  the  first  problem  above.  The  technique 
determines  whether  particles  of  low  weight  are  allowed  to  continue  the  random 
walk.  Implementation  of  Russian  roulette  requires  a  minimum  stipulated  weight 
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Figure  I.  Random  Walk  Flow-Chart 


value  MORbE  uses  to  determine  which  particles  undergo  the  technique.  Addition 
ally.  MORSE  requires  an  average  weight  value  to  assign  to  particles  surviving 
Russian  roulette.  The  chance  of  survival  is  the  ratio  of  the  lower  weight  value  to 
the  average  weight  value  (I6:29d). 

Splitting  is  a  useful  technique  to  increase  the  number  of  particles  in  regions 
more  than  a  few  mean  free  paths  from  the  surface  of  the  geometry.  Within  a 
user-specified  region.  MORSE  splits  a  particle  if  the  weight  is  greater  than  a 


specified  minimum  value.  The  first  step  is  multiplication  of  the  particle  weight  by 
E^/E,  to  account  for  the  fraction  lost  to  absorption.  Next,  MORSE  divides  the 
weight  by  two  to  account  for  the  two  particles  produced.  The  splitting  continues 
until  the  particle  weights  fall  below  the  minimum  value  mentioned  above.  This 
method  produces  many  more  statistical  particles  possibly  resulting  in  a  reduction 
of  sample  variance  (4:16). 

In  order  to  use  MORSE  for  a  specific  problem,  problem-specific  input  data 
and  problem  dependent  user-modified  subroutines  must  be  produced.  The  input 
data  includes  parameters  controlling  Monte-Carlo  random  walk,  combinatorial 
geometry  description,  cross  section  information,  and  user-requested  calculations. 
The  user-modified  subroutines  (specific  to  this  thesis)  include  routines  to  generate 
source  neutrons,  assign  cross  section  media  regions,  determine  dose  contributions 
from  single  neutron  events,  and  read  in  additional  input  data. 

Human  Phantom 

The  human  phantom  created  in  this  project  is  5  feet,  9.3  inches  (176  cm)  tall 
with  a  mass  of  151  ±  2  pounds  (69  ±  1  kg).  The  phantom  is  modeled  from  the 
ICRP  reference  man  (10),  Gray’s  anatomy(9),  the  Snyder  phantom  model 
(8:1473-1478),  and  discussions  with  a  medical  doctor  (14).  The  reference  man  is  5 
feet,  9  inches  (175  cm)  tall  and  weighs  154  pounds  (70  kg).  The  mass  of  the 
phantom  model  differs  from  that  of  the  reference  man  primarily  because  the  phan¬ 
tom  is  somewhat  more  simplistic  than  the  reference  man.  Figure  2  shows  a 
three-dimensional  representation  of  the  model  used  in  this  study.  The  figure  on 
the  left  represents  the  skin  of  the  phantom.  The  figure  on  the  right  shows  the 
phantom  bones. 


Figure  2.  Three-Dimensional  Representation  of  Human  Phantom 


The  coordinate  axis  of  the  phantom  originates  at  the  center  of  the  base  of  the 
right  elliptical  cylinder  defining  the  torso.  The  positive  x-direction  is  toward  the 
left-hand  side  of  the  body,  the  positive  y-direction  is  toward  the  back  of  the  body, 
and  the  positive  z-direction  is  toward  the  head.  Sixtv-one  geometric  shape  entries 
are  combined  to  produce  the  phantom  with  its  five  organs:  «kin,  bone,  brain, 
gastro-intestinai  tract,  and  muscle  (includes  the  remaining  tissue  of  the  body). 
Additional  organs  are  not  modeled  since  this  thesis  is  (primarily  interested  in) 


determining  incapacitating  doses  to  the  body.  The  actual  data  lines  that  produce 


the  phantom  are  included  in  the  input  data  file  listed  in  Appendix  A.  Figures  3 


and  4  are  computer  generated  phantom  representations  to  verify  correctness  of  the 


input.  The  first  representation  (Figure  3)  is  a  cross  sectional  front  view  of  the 


phantom  cut  down  the  center  plane  (y=0).  The  darker  points  represent  bone,  the 


smaller  points  represent  muscle,  the  dashes  are  the  gastro-intestinal  tract,  and  the 


dark  boxes  represent  the  brain.  The  relatively  small  thickness  of  the  skin  does 


not  allow  for  a  computer  representation  on  this  scale.  The  second  (Figure  4)  is  a 


cross  sectional  side  view  cut  down  the  center  plane  (x=0).  The  representations  of 


Figure  4  are  the  same  as  for  Figure  3.  Additionally,  the  skin  in  representable  on 


this  size  view  (shown  as  "+"  signs).  The  representation  in  Figure  4  is  somewhat 
"squashed"  down  but  accurately  depicts  the  relative  position  of  each  organ.  Both 


representations  are  produced  using  an  RSIC  code  written  specifically  to  create 


such  images  (13:Sect  3). 


The  input  data  file  in  Appendix  A  contains  the  exact  descriptions  of  the  62 


geometric  shapes  required  to  produce  the  phantom.  The  individual  depictions  of 


each  phantom  component  are  included  in  the  thesis  notebook  (1).  The  phantom 


model  consists  of  four  main  regions:  head,  torso,  legs,  and  arms. 


Head.  The  head  is  a  sphere  of  outer  radius  9.0  cm,  with  a  right  circular 


cylinder  of  radius  7.5  cm  attached  to  form  the  jaw.  Another  right  circular  cylinder 


of  radius  5.2  cm  attaches  from  below  the  sphere  to  form  the  neck.  The  skin  covers 


the  head  with  a  thickness  of  0.2  cm  (same  thickness  over  the  entire  phantom). 


The  skull  is  0.5  cm  in  thickness  adjacent  to  the  skin  of  the  head.  The  spinal  cord 


begins  at  the  base  of  the  skull  and  extends  to  the  torso  as  described  in  the  next 


paragraph.  The  brain  is  approximately  two-thirds  of  a  sphere  of  radius  8.3  cm 


within  the  spherical  portion  of  the  skull.  The  remaining  space  within  the  head  is 


muscle. 


V? 


Torso.  The  torso  is  a  right  elliptical  cylinder  70  cm  tall  with  a  minor  radius 
of  10.5  cm  and  a  major  radius  of  15.5  cm.  Beneath  the  skin  lies  the  ribcage,  spinal 
cord,  pelvis,  and  gastro-intestinal  tract.  The  ribcage  is  a  right  elliptical  cylindri¬ 
cal  shell  cut  by  cross  sectional  planes  to  produce  11  ribs  0.5  cm  deep  and  1.7  cm 
high  (see  Figure  2).  The  lower  five  ribs  are  cut  by  wedges  to  produce  an  abdomi¬ 
nal  opening.  Additionally,  a  0.5  cm  thick  sternum  is  over  the  top  six  ribs. 

The  spinal  cord  is  a  2.0  cm  thick  right  circular  cylinder  originating  beneath 
the  skull  and  angling  toward  the  middle  portion  of  the  ribcage.  From  there,  the 
spinal  cord  follows  the  ribcage  straight  down  to  the  pelvis.  The  pelvis  is  a  right 
elliptical  cylindrical  shell  with  a  portion  cut  off  in  front  (see  Figure  2).  The  pelvis 
originates  at  the  base  of  the  torso  and  is  19.6  cm  high.  The  gastro-intestinal  tract 
is  a  right  elliptical  cylinder  originating  inside  and  approximately  half  way  up  the 
pelvis.  The  tract  extends  an  additional  5.2  cm  above  the  pelvis.  All  additional 
space  within  the  torso  is  muscle. 

Legs.  The  legs  are  divided  into  upper  and  lower  legs,  both  truncated  right 
cones.  The  upper  and  lower  radii  of  the  entire  leg  are  8.5  cm  and  3.5  cm,  respec¬ 
tively.  Within  each  leg  are  the  upper  and  lower  leg  bones.  The  upper  leg  bones 
are  right  circular  cylinders  of  radius  2.5  cm.  The  bones  originate  at  the  outer  base 
of  the  pelvis  and  extend  to  the  center  location  between  the  upper  and  lower  legs. 
The  lower  leg  bones  (tibia  and  fibula)  are  modeled  as  a  single  right  circular 
cylinder  within  each  leg. 

Arms.  The  arms  are  also  divided  into  upper  and  lower  sections.  The  radii  of 
the  upper  and  lower  portions  of  the  arms  are  4.0  and  2.0  cm,  respectively.  The 
arm  bones  are  right  circular  cylinders  of  1.5  cm  (upper  arm)  and  1.25  cm  (lower 
arm  -  ulna  and  radius  combined)  through  the  center  of  each  arm. 

All  of  the  organs,  with  exception  of  the  bone,  are  approximated  as  homogene¬ 
ous  tissue  with  a  density  of  0.987  g  / cm3  (10:290-324).  Although  this  may  seem  to 


be  a  crude  approximation,  actual  organ  densities  are  very  difficult,  if  not  impossi¬ 
ble,  to  determine.  Throughout  the  day  organs  vary  in  composition  because  of  typ¬ 
ical  body  functions.  However,  the  average  daily  mass  to  volume  ratio  of  most 
organs  is  almost  exactly  that  of  tissue  (14).  The  bone  is  modeled  with  a  density 
of  1.49  g /cm3.  The  constituents  of  the  phantom  tissue  and  bone  are  listed  in 
Table  I.  Table  II  presents  the  phantom  organ  volumes  and  their  uncertainty  in 
measurement.  The  uncertainties  result  since  the  more  complex  geometric  combi¬ 
nations  produce  non-trivial  volume  equations. 

Cross  Section  Mixing 

The  MORSE  computer  package  includes  a  choice  for  mixing  cross  sections: 
mix  during  each  problem  analysis  or  pre-mix  prior  to  problem  analysis.  The 
phantom  problem  requires  the  latter  because  of  buffer  storage  problems  resulting 
from  the  large  storage  requirements  of  the  geometry  and  particle  histories. 

Three  media  are  required  for  this  problem:  soft  tissue,  bone,  and  air.  Table  I 
contains  the  elemental  constituents  of  soft  tissue  and  bone.  The  air  is  a  mixture 
of  nitrogen  (80.0  percent  by  weight)  and  oxygen  (20.0  percent  by  weight).  The 
individual  elemental  cross  sections  are  from  the  RSIC  DLC-37  data  library  con¬ 
taining  37  neutron  and  21  gamma-ray  energy  group  cross  sections  (2).  The 
library  is  an  ANISN  format  library  with  the  following  cross  sections:  absorption, 
fission,  total,  within  group,  and  group-to-group.  The  scattering  coefficients  are 
represented  by  a  P3  Legendre  polynomial  expansion  of  the  differential  scattering 
cross  sections. 
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Figure  3.  MORSE  cross  sectional  Representation  of  Phantom  from  Front 
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Figure  4.  MORSE  cross  sectional  Representation  of  Phantom  from  Side 
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Evolution  of  Final  MORSE  Version 
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The  final  version  of  MORSE  determines  normalized  absorbed  radiation  dose 
to  a  human  phantom  for  any  of  the  five  source  descriptions:  isotropic  and  mono- 
directional  from  right  side,  above,  front,  and  back.  By  simply  modifying  the 
input  data  file,  a  wide  range  of  problems  can  be  solved  without  any  modifications 
to  MORSE.  Appendix  A  contains  an  example  input  data  file.  The  remainder  of 
this  section  covers  the  10  stages  of  MORSE  evolution  resulting  in  the  final 
MORSE  version  of  this  thesis. 

Familiarization.  The  MORSE  code  package  at  the  Air  Force  Institute  of 
Technology  (AFIT)  is  a  version  modified  to  run  on  a  Digital  Electronics  Corpora¬ 
tion  VAX/VMS  system.  Subsequently,  none  of  the  many  FORTRAN  files  require 
modification  if  the  code  is  run  on  the  AFIT  CSC  computer  system.  However,  as  a 
check  to  ensure  program  integrity,  the  code  includes  sample  problems  complete 
with  results.  The  results  of  the  first  sample  problem  (13:  Sect  11,2-1)  run  at  AFIT 
did  indeed  match  those  forwarded  by  RSIC. 

Initial  Test.  This  stage  covers  the  initial  modification  of  MORSE  subroutines 
and  a  comparison  with  diffusion  equation  results.  The  problem  is  a  sphere  of 
water  with  a  15  cm  radius.  The  sphere  receives  an  isotropic  neutron  fluence  of 
1  n /cm  2  at  the  sphere  surface  (directed  over  27T  steradians  inward).  The  cross 
sections  are  for  thermal  energy  neutrons  in  water  but  are  used  over  a  range  of  1 
Mev  to  0.025  ev  in  this  1-group  problem.  The  initial  test  includes  rewriting  the 
fluence  estimator  subroutines  to  include  a  track  length  estimator  for  spherical 
geometry.  The  track  length  estimator  equation  is: 


where  </>  is  the  normalized  fluence  in  n /cm2/ source  n.  /,  is  t he  distance  traveled 
(tracklength)  by  particle  i  in  the  sphere  (cm),  N  is  the  number  of  particles  tracked 
by  MORSE,  and  V  is  the  volume  of  the  sphere  (cm3).  MORSE  models  the 
incident  neutron  fluence  on  the  sphere  surface  using  random  number  generators  to 
direct  the  neutrons  isotropically  inward.  MORSE  also  determines  whether  or  a 
not  a  neutron  is  within  the  sphere  volume  using  a  simple  test:  is  the  distance 
between  the  neutron  and  the  center  of  the  sphere  greater  than  the  radius  of  the 
sphere?  If  the  distance  is  larger,  the  neutron’s  tracklength  is  not  included  in  the 
fluence  estimate.  A  sample  population  of  5000  source  neutrons  produces  a  fluence 
estimate  in  the  sphere  volume  of  4.0  X  10-4n  / cm2 / source  n  with  a  standard 
deviation  of  5.0%. 

The  solution  to  the  diffusion  equation  (1:29)  for  the  water  sphere  problem  is: 

0.71  sinh(Br) 
r 

where  <z5(r )  is  the  fluence  at  radius  r  (n / cm2),  and  B  is  the  material  buckling 
(c/n-1).  The  diffusion  equation  solution  is  for  a  non-multiplying  medium  using 
one-group  thermal  neutron  cross  sections.  Integrating  eq.  (2)  over  the  sphere’s 
radius  and  dividing  by  the  volume  produces  a  volume  average  fluence  over  the  15 
cm  radius  sphere  of  0.37  n /cm2  (1:32).  In  order  to  compare  solutions  between 
MORSE  and  the  diffusion  equation,  the  MORSE  results  are  multiplied  by  the 
total  number  of  source  neutrons  from  a  normalized  fluence  of  1  n /cm2  on  the 
sphere  surface  {Ai^lbcm)2  * \n  / cm2  —  2830  source  neutrons)  to  produce  1.1 
n /cm2.  However,  the  MORSE  results  are  still  a  factor  of  two  larger  since  the 
source  model  was  only  over  2n  and  not  47r  steradians.  The  final  MORSE  result  is 
then  0.55  n /cm2  ±  0.03.  The  MORSE  and  diffusion  equation  results  are  within 
an  order  of  magnitude  of  each  other. 


(2) 


Multiple  Volume  Test.  The  next  MORSE  modification  solves  the  previous 
problem  with  a  slight  variation:  the  source  is  on  a  spherical  surface  with  a  radius 
of  35  cm,  concentric  to  the  original  sphere.  This  modification  is  a  benchmark  for 
future  geometry  modifications.  The  MORSE  analysis  of  the  2000  neutrons 
directed  isotropically  inward  results  in  a  normalized  fluence  in  the  inner  sphere  of 
4.4  X  10 ~sn/cm~/ source  n  with  a  standard  deviation  of  7.0%. 

Non-Spherical  Geometry.  The  modifications  prior  to  this  stage  are  strictly 
for  concentric  sphere  problems.  The  modification  of  this  stage  is  to  determine 
fluence  in  any  three-dimensional  shape  using  geometry  independent  tests.  The  new 
tests  determine  the  current  region  of  the  neutron  location.  A  region  is  a  user- 
defined  volume  described  in  the  combinatorial  geometry  section  of  the  input  data 
file.  Using  the  same  problem  as  in  the  multiple  volume  test  but  with  the  new 
subroutine  modification,  MORSE  again  produces  a  normalized  fluence  of 
4.4  X  I0~s n  /  cm  2 / source  n.  The  result  validates  the  geometry  independent 
modification  in  the  spherical  problem. 

Non-Spherical  Test.  This  stage  determines  the  validity  of  the  previous 
modification  on  non-spherical  geometry.  The  problem  consists  of  two  cubes,  full  of 
vacuum,  as  shown  in  Figure  5.  Neutrons  are  incident  on  the  left-hand  surface  of 
cube  1  in  the  +x  direction.  The  only  modification  to  MORSE  is  a  statement  after 
each  source  neutron  is  created  to  print  the  starting  location.  Knowing  the  starting 
location  determines  whether  a  neutron  travels  only  through  cube  1  or  both  cube  1 
and  cube  2.  Since  the  cubes  are  full  of  vacuum,  the  neutrons  do  not  collide  and 
change  direction.  The  code  produces  normalized  fluences  of  0.250n /cm2 /source  n 
and  0.0875n /cm 2 / source  n  for  cube  1  and  cube  2,  respectively,  from  20  starting 
source  neutrons.  Hand  calculation  (1:39)  verifies  these  results  since  all  20  must 
pass  through  cube  1  while  only  7  pass  through  cube  2. 
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Figure  5.  Geometry  for  Non-Spherical  Test. 

Multiple- Group  Upgrade.  The  next  stage  is  the  conversion  from  a  single 
group  water  problem  to  a  37  group  tissue  problem  using  the  geometry  defined  in 
Stage  2  (source  on  surface  of  15  cm  radius  sphere).  An  analysis  of  2000  neutrons 
produces  a  normalized  fluence  in  the  sphere  of  1.52  X  1 0-3n /cm2/ source  n  with  a 
standard  deviation  of  o.3%.  This  fluence  is  approximately  three  times  larger  than 
the  fluence  produced  using  one-group  cross  sections;  the  difference  is  attributable 
to  the  ability  of  each  neutron  to  scatter  from  group  to  group,  with  varying  cross 
sections,  between  1  Mev  and  0.025  ev  (group  21  to  group  37). 

Phantom  Geometry  Addition.  The  seventh  stage  is  the  first  to  use  the  phan¬ 
tom  geometry.  The  main  purpose  of  this  stage  is  to  determine  fluence  sensitivity 
to  change  in  organ  volume.  Table  II  shows  three  of  the  five  organs  (skin.  bone, 
and  muscle)  have  an  uncertainty  in  volume  measurement.  To  determine  the 


f 


sensitivity  of  the  estimated  fiuence  to  the  volume  uncertainty,  five  analyses  are 


run  through  MORSE.  In  each  case,  only  the  organ  volumes  are  changed  in  the 
input  data  file  since  the  dose  effect  from  individual  tracklengths  is  divided  by  the 
volume  (eq,  1).  Table  III  shows  the  results  of  the  five  1000-neutron  analyses. 

Table  III 

Fiuence  Sensitivity  to  Organ  Volume  Uncertainty 


Normalized  Fiuence  [n/cm* /source  nl 

Skin 

Bone 

Muscle 

9.46E-6 

1.16E-5 

9.94E-6 

9.30E-6 

1.16E-5 

9.94E-6 

9.46E-6 

1.12E-5 

9.94E-6 

9.46E-6 

1.16E-5 

9.85E-6 

9.30E-6 

1.12E-5 

9.85E-6 

The  first  analysis  uses  the  mean  phantom  organ  volumes  shown  in  Table  II. 
The  second  analysis  uses  the  upper  limit  of  the  skin  volume  and  leaves  the  bone 
and  muscle  unchanged.  The  third  analysis  uses  the  upper  limit  of  the  bone  volume 
and  mean  volumes  of  skin  and  muscle.  The  fourth  analysis  uses  the  upper  limit  of 
the  muscle  and  mean  volumes  of  the  skin  and  bone.  The  fifth  analysis  uses  the 
upper  limit  of  each  organ  volume.  The  results  show  that  the  skin's  volume  uncer¬ 
tainty  varies  the  fiuence  by  no  more  than  1.7%;  the  bone  volume  uncertainty 
varies  the  fiuence  by  no  more  than  3 A%;  and  the  muscle  volume  uncertainty 
varies  the  fiuence  by  no  more  than  0.9%.  The  fifth  analysis  shows  that  whether 
the  organ  volumes  are  varied  singly  or  all  at  once,  their  net  effect  is  the  same:  a 
variation  of  a  specific  organ  volume  only  affects  that  organ's  fiuence.  The  remain¬ 
ing  two  organs  have  no  volume  uncertainties  and  do  not  show  a  variation  in 


fiuence  results  from  the  above  analyses. 


Absorbed  Dose  Calculation.  The  next  stage  is  the  conversion  of  MORSE  out¬ 
put  from  normalized  fluence  to  normalized  absorbed  dose.  The  conversion  method 
recommended  by  the  International  Commission  on  Radiation  Units  and  Measure¬ 
ments  (15:95)  is  use  of  KERMA  factors.  KERMA  stands  for  Kinetic  Energy 
Released  in  MAterial  and  is  the  sum  of  initial  kinetic  energies,  of  charged  particles 
resulting  from  neutron  interaction  in  a  volume,  per  mass  of  material  in  the 
volume.  KERMA  factors  are  actually  a  KERMA  to  unit  fluence  ratio.  The 
KERMA  factors  used  in  this  project  are  from  the  RSIC  cross  section  library  men¬ 
tioned  above  in  the  cross  section  mixing  section.  Individual  elemental  KERMA  for 
the  11  constituent  elements  of  tissue  and  bone  are  mixed  to  produce  tissue  and 
bone  KERMA  factors  for  all  37  energy  groups.  The  tissue  and  bone  KERMA  fac¬ 
tors  are  listed  in  Table  IV  (energy  ranges  of  the  groups  are  listed  in  Table  V). 

The  normalized  dose  for  organ  m  is  then  determined  as: 

Dtm  =  (j>™  *  KF™  [RAD— {tissue  or  bone)/{source  n)].  (3) 

where  <$>/”  is  the  MORSE  estimate  of  normalized  fluence  (Eq.  1)  of  energy  group  i 
and  KF,m  is  the  KERMA  factor  of  the  ith  energy  group. 

A  total  absorbed  dose  is  calculated  by  multiplying  the  MORSE  output  of 
normalized  dose  with  the  total  incident  fluence: 

Dm  =  Dm  *  4>inc  (4) 

where  Dm  is  the  total  absorbed  dose  of  organ  m  and  <f>inc  is  the  incident  neutron 
fluence. 

Sample  Analysis.  The  ninth  stage  of  the  MORSE  evolution  is  analysis  of  the 
phantom  using  a  fission  weapon  energy  spectrum  (2:11)  with  neutrons  incident 
five  different  ways: 

1)  isotropic  from  a  sphere  surrounding  the  phantom  (radius  =  100  cm, 
center  at  0,0,10), 


2)  mono-directional  in  the  +x  direction  (incident  on  the  right  side  of  the 
phantom  from  a  source  plane  defined  by  x  =  —30,  —13  <  y  <13, 

-80  <  2  <  96), 

3)  mono- directional  in  the  -z  direction  (incident  from  above  the  phantom 
from  a  source  plane  defined  by  —30  <  x  <  30,  —13  <  y  <  13,  2  =  96), 

4)  mono-directional  in  the  +y  direction  (incident  on  the  front  of  the  phan¬ 
tom  from  a  source  plane  defined  by  —30  <  x  <  30,  y  =  — 13, 

-80  <  2  <  96), 

5)  mono-directional  in  the  -y  direction  (incident  on  the  back  side  of  the 
phantom  from  a  source  plane  defined  by  —30  <  x  <  30,  y  =  13, 

-80  <  2  <96. 

The  results  are  presented  in  Chapter  III,  Results.  Table  V  shows  the  energy  distri¬ 
bution  of  the  weapon  spectrum. 

Final  MORSE  Version.  The  final  stage  is  the  creation  of  the  user-ready  ver¬ 
sion  of  MORSE.  Prior  to  this  stage,  using  different  source  directions  required 
modification  to  the  MORSE  FORTRAN  subroutines.  Instead,  MORSE  is  modified 
to  read  special  input  data  describing  the  source  direction  required.  Additionally, 
the  surface  area  of  the  respective  source  type  is  multiplied  with  the  dose  per 
source  neutron  to  produce  dose  per  unit  fluence  (normalized  dose)  in  units  of 
GRAY-(tissue  or  bone) / (n  / cm2).  Appendix  B  contains  the  final  modified  MORSE 
subroutines.  The  results  of  this  stage  are  presented  in  Chapter  III  of  this  thesis. 


Table  1\  . 


Tissue  and  Bone  IvERMA  Factors 


KERMA  Factors 

Tissue 

Bone 

[RAD-Tiss /(n  /  cm2)] 

[RAD-Bone/U  .  cm2)] 

0.723zio^* 

0.551  z  io-08 

0.686z  to-0® 

0.520z  io-® 

0.664x  io-® 

0.502z  io^* 

0.653x  10“® 

0.493x10“® 

0.638x  io-® 

0.482xio“® 

0.619z  nr08 

0.468ziO“® 

0.612z  10“® 

0.461x10“® 

0.586ziO“® 

0.440xio“® 

0.569x10“® 

0.427x10“® 

0.539ziO“® 

0.399ziO“® 

0.528z  10"® 

0.391Z10"08 

0.504ziO“® 

0.366x10“® 

0.472xio*<* 

0.34  lxio^* 

0.448x10“® 

0.320x10“°* 

0.441x10“® 

0.316x10“® 

0.419x10“® 

0.300z  io-08 

0.356ziO“® 

0.249xio“® 

0.326z  io-® 

0.226z io-® 

0.320xiO“® 

0.221z  i0“® 

0.276ziO“® 

0.191xio“® 

0.210xiO“® 

0.145ziO“® 

0.130z  io-® 

0.894x10“® 

0.801xiO“® 

0.547ziO“® 

0.559x10^° 

0.381xiO“® 

0.315ziO“® 

0.214xiO“® 

0.215x10“® 

0.146x10“® 

0.147ziO“® 

O.lOOxio-® 

0.627xiO“‘° 

0.428xiO“‘° 

0.221xio“‘° 

0.152x10“'° 

0.919zio-‘l 

0.645x10“" 

0.310zio-‘‘ 

0.241xiO““ 

0.103zio“" 

0.120xio“" 

0.919ziO“12 

0.141x10“" 

0.134z  io-" 

0.222iiO“" 

0.225zio-“ 

0.375x10“" 

0.369zio-“ 

0.610xiO“" 

0.151x10“'° 

0.247x10“'° 

Table  V. 
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Fission  Weapon  Neutron  Energy  Spectrum 
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1.96a:  10+07 
1.69xlO+07 
1.49xlO+07 
1.42xlO+07 
1.38xlO+07 
1.28x  10+07 
1.22xlO+07 
l.llxlO+07 
1.00xl0+°7 
9.05a:  10+06 
8.19xlO+06 
7.41xlO+06 
6.38xlO+06 
4.97xlO+06 
4.72xlO+06 
4.O7xlO+06 
3.01  xlO+06 
2.39xlO+06 
2.31xlO+06 
1.83xlO+06 
l.llxlO+06 
5.55x  10+05 
1.58x  10+05 
l.llxlO+05 
5.25x  10+°4 
2.48xlO+04 
2.19x  10+04 
1.03x  10+04 
3.35x  10+03 
1.23xlO+03 
5.83xlO+02 
1.01xl0+°2 
2.90xl0+f)1 
1.07xl0+01 
3.06x  10+°° 
1.13x  10+°° 
4.14x  10~Q1 


F  raction 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.004 

0.004 

0.005 

0.007 

0.018 

0.003 

0.009 

0.055 

0.032 

0.011 

0.097 

0.147 

0.216 

0.150 

0.019 

0.121 

0.057 

0.006 

0.024 

0.014 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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III.  Results 


Overview 

The  user-ready  version  of  MORSE  produced  from  the  last  stage  of  the  code 
evolution  is  the  version  used  for  analysis  in  this  chapter.  The  results  are  by  no 
means  for  an  actual  problem  but  instead  show  how  MORSE  is  used,  in  its 
modified  form,  to  determine  absorbed  neutron  dose  to  organs  of  a  phantom. 

This  chapter  contains  four  results  sections;  analyses  using  the  user-ready  ver¬ 
sion  of  MORSE,  a  comparison  of  MORSE  results  to  results  from  another  study, 
variance  and  processing  time  reduction  analyses,  and  a  sample  calculation  using 
MORSE  output  to  determine  estimated  dose  to  a  soldier.  In  all  cases,  the  analyses 
consist  of  neutrons  incident  on  the  phantom  in  each  of  the  five  source 
configurations  described  in  the  Sample  Analysis  Stage  of  Chapter  II.  The  energy 
spectrum  is  the  fission  weapon  spectrum  shown  in  Table  V.  Because  of  the  buffer 
storage  problem  described  in  Chapter  TI,  MORSE  is  capable  of  analyzing  only 
1200  particles  per  batch  in  problems  with  no  splitting  and  1000  particles  per 
batch  in  problems  with  splitting. 

Test  Case  Results. 

The  first  analysis  (referred  to  as  "test  case")  uses  a  sample  population  of 
24,000  neutrons  (20  batches  of  1200  neutrons  per  batch).  The  results  are 
presented  in  Table  VI.  In  each  box  of  Table  VI,  the  top  number  is  the  normalized 
(unit)  dose  and  the  bottom  number  is  the  fractional  standard  deviation.  The 
table  shows  the  form  of  results  attainable  for  a  given  neutron  energy  spectrum. 
Multiplication  of  any  of  the  organ-source  direction  results  with  an  incident  fluence 
provides  a  total  absorbed  dose  to  the  organ  (eq.  4).  A  solution  set  such  as  Table 
VI  eliminates  the  need  for  any  future  computer  analysis  for  a  given  energy  spec¬ 
trum.  Note  the  large  central  processing  unit  (CPU)  times  required  for  the  mono- 


:v>.v 


directional  sources  compared  with  the  isotropic  source.  These  larger  CPI.  times 
occur  since  many  more  source  neutrons  interact  with  the  phantom  than  in  the  iso¬ 
tropic  case  (see  source  plane  descriptions  in  the  Sample  Analysis  Stage  of  Chapter 


Table  VI. 

Normalized  Dose  Results  -  Test  Case 


Organ  # 


Bone 


Muscle 


Brain 


GI  Tract 


CPC  (min) 


Norm.  Dose  fGrav/fn  /cm "II  from  given  source 


ISOTROPIC 


6.99:r  10-12 
0.05 


3.662  10“12 
0.07 


4.272 10“12 
0.06 


4.912 10"12 
0.15 


1.222 10-12 

0.39 


12.70 


Rt.  Side 


1.09210-11 

0.01 


5.402 10-12 

0.02 


6.432 10-12 

0.01 


8.482 10"12 
0.05 


2.662 10-12 
0.15 


84.76 


Above 


3.53210-12 

0.02 


1.102 10-12 
0.02 


1.222 10“ 12 
0.03 


1.122 10-11 
0.02 


4.I62 10~15 

0.87 


46.43 


Front 


1.47210-11 

0.02 


7.36210-12 

0.03 


1.172l0~u 

0.02 


5.422 10“ 12 
0.08 


6.662IO'12 

0.08 


80.90 


Back 


9.41210-12 

0.02 


1 . 1 1 2 10~ 1 1 
0.01 


1.172 10-u 
0.08 


6.762 10 

0. 


Note:  the  numbers  below  normalized  dose  are  fractional  standard  deviation. 


Comparison  of  Results. 


The  vast  majority  of  open  literature  on  neutron  doses  addresses  results  in 
terms  of  transmission  factors.  A  transmission  factor  is  the  ratio  of  KERMA  in  an 
organ  to  soft-tissue  KERMA  in  free  air  (no  organ  present).  The  organ  KERMA 
for  the  phantom  problem  is  simply  the  MORSE  output,  since  the  output  is  a  nor¬ 
malized  fluence:  KERMA  =  GRAY /(n  /  cm2)  *  1  n/cm2  —  GRAY.  The 
KERMA  in  free  air  (more  commonly  referred  to  as  "free-in-air  KERMA")  for  the 
phantom  problem  is  determined  (1:58)  as: 

K'fia  =  KF tissue  *  f'  (5) 

which  leads  to 

K/l.  -  £*},.  (6) 

1-1 

where  KFlissue  is  the  tissue  KERMA  factor  of  the  ith  energy  group  and  / 1  is  the 
source  fission  energy  fraction  of  the  ith  group.  Table  VII  fists  the  free-in-air 
KERMA  values  for  all  37  neutron  energy  groups.  The  total  free-in-air  tissue 
KERMA  is  2.1  X  10-9  RAD-tissue. 

The  only  two  organs  applicable  for  comparison  with  results  from  other  stu¬ 
dies  are  the  brain  and  skeleton.  This  is  primarily  because  the  remaining  organs  of 
the  human  phantom  studied  here  are  not  commonly  analyzed  by  others  (this 
study  is  primarily  interested  in  immediate  incapacitation  to  organs).  Table  VIII 
shows  a  comparison  of  the  MORSE  bone  and  brain  phantom  transmission  factors 
to  transmission  factors  of  Japanese  atomic  bomb  survivors  (6:617).  The  Japanese 
survivor  data  is  from  a  joint  Science  Applications  International  Corporation 
(SAIC)  and  RSIC  data  base. 
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Table  MI 


30 

31 

32 

33 

34 

35 

36 

37 


Total 


Free-in-Air  KERALA  \  alues 


F  Faction 

IvER\LA  Factor 
RAD-Tis.  ( n  / cm  "11 

0.0000 

0.723xl0"8 

0.0000 

0.686x  10“8 

0.0000 

0.664x  10“8 

0.0000 

0.653xl0"8 

0.0000 

0.638xl0"8 

0.0000 

0.619xl0~8 

0.0000 

0.612xl0“8 

0.0000 

0.586xl0-8 

0.0038 

0.569xl0"8 

0.0035 

0.539xl0"8 

0.0054 

0.528x  10~8 

0.0074 

0.504xl0-8 

0.0184 

0.472xl0-8 

0.0033 

0.448x  10-8 

0.0085 

0.441x  10-8 

,n 

II  l|ll.(  ID  X 

IUM2  l 

l).:’.56x  It)-3 

0.0106 

0.326xl0“8 

0.0972 

0.320xl0“8 

0.1470 

0.276x  10-8 

0.2160 

0.2  lOx  10-8 

0.1500 

0.1 30x  10-8 

0.0193 

O.SOlxlO-9 

0.1210 

0.559xl0-9 

0.0573 

0.315x  10-9 

0.0060 

0.215x  10~9 

0.0240 

0.1 47x  10-9 

0.0144 

0.627  x  10“ 

0.0000 

0.221xl0-1° 

0.0000 

0.919x  10-1 1 

0.0000 

0.310xl0_u 

0.0000 

0.103x10" 11 

0.0000 

0.919X10"1- 

0.0000 

0.134X10"11 

0.0000 

0.225x10“ 11 

0.0000 

0.369x10" 11 

0.0000 

0.151x10“ 10 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.219x10 
0.189x10 
0.285x10 
0.370x10 
0.869x10 
0.146x10 
0.374x10 
it  -':t |  i  in 
o.i  15./  I() 
0.346x10 
0.311x10 
0.406x10 
0.453x10 
0.195x10 
0.155x10 
0.676x10 
0.180x10 
0.129x10 
0.353x10 
0.903x10 


1.1 00x10' 


Table  Mil 


IvERMA  Transmission  Factors:  Phantom  versus  Japanese  Survivors 


MORSE 

SAIC/ORNL 

Bone 

0.37 

0.36 

(0.03) 

(0.05) 

Brain 

0.31 

0.37 

(0.12) 

10.051 

The  numbers  in  parentheses  are  fractional  standard  deviations. 

The  transmission  factor  comparison  is  very  close  for  the  bone  results.  The 
brain  comparison  is  not  as  good  and  therefore  prompted  another  look  at  the 
modeling  of  the  head  (from  this  study).  The  size  of  the  brain  from  the 
SAIC/ORNL  model  is  not  known  for  comparison,  but  the  phantom  brain  from 
this  study  may  be  approximately  100  cm3  too  large  if  the  Snyder  model  is  accu¬ 
rate  (8:1476).  As  a  check,  a  reduction  in  the  size  of  the  phantom  brain  by  100cm3 
produces  a  transmission  factor  of  0.32.  Any  further  discrepancies  in  the  transmis¬ 
sion  factors  cannot  be  fully  realized  without  knowing  the  details  of  the 
SAIC/ORNL  model.  One  can  assume  the  skull  thickness  (unlfcuui  in  this  study) 
attenuates  too  many  neutrons  resulting  in  a  lower  transmission  factor. 

Time/Variance  Reduction  Results. 

This  third  section  of  results  contains  four  attempts  to  reduce  sample  variance 
and/or  CPU  time.  The  first  series  uses  splitting,  Russian  roulette,  and  time  kill 
for  24.000  source  neutrons  (24  batches  of  1000).  Time  kill  is  a  function  that 
creates  a  neutron  "death"  for  neutrons  older  than  10  microseconds.  Ten 
microseconds  is  chosen  to  ensure  the  neutrons  reach  the  lowest  energy  group 
before  being  eliminated.  The  results  are  presented  in  Table  IX. 
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The  second  series  uses  splitting  and  time  kill,  again  for  2-1  batches  of  1000 
neutrons.  The  results  are  presented  in  Table  X. 

The  results  of  the  third  series  are  presented  in  Table  XI.  The  series  uses  only 
time  kill  for  20  batches  of  1200  particles. 

The  fourth  and  final  series  uses  only  time  kill  but  runs  40  batches  of  1200 
neutrons.  The  results  are  presented  in  Table  XII. 


Table  IX 

Normalized  Dose  Results  -  Series  #1 


Organ  # 

Norm.  Dose  [Grav/(n /cm 2)]  from  eiven  sonrcp 

ISOTROPIC 

Rt.  Side 

Above 

Front 

Back 

Skin 

6.93x  10“12 
0.05 

1.06x  10"11 
0.02 

3.48xl0"12 

0.02 

1.16x10"" 

0.01 

1.54x10" 11 
0.02 

Bone 

3.39xl0“12 

0.06 

5.13xl0“12 

0.02 

1.19xlO"12 

0.03 

7.34xl0"12 

0.02 

9.92x  10~12 
0.02 

Muscle 

4.26i  10-12 
0.05 

5.13x  10"12 
0.02 

1.25x  10~12 
0.02 

1.12x10"" 

0.02 

1.10x10"" 

0.02 

Brain 

5. 45x10" 12 
0.17 

9.58x  10" 12 
0.05 

1.07x10"" 

0.03 

o.64x  10~12 
0.13 

9.95x10" 12 
0.10 

GI  Tract 

1.42x10" 12 
0.31 

2.35x10" 12 
0.11 

2.73xl0~14 

0.49 

7.59x  10"12 
0.08 

6.03x10" 12 
0.07 

CPU  (min’!  I 

1 1 .55 

63.69 

12.57 

71.96 

69.07 

Table  \ 


Normalized  Dose  Results  -  series 


Organ  # 


Bone 


Muscle 


Brain 


GI  Tract 


CPU  (min 


Norm.  Dose  fGrav  i  n  7 cm  ■'ll  from  givpn  source 


ISOTROPIC  Rt.  Side 


6.60x  10“ 12 
0.045 


3.60x  10“12 
0.081 


3.96xl0“12 

0.044 


5.07x10  12 
0.234 


1.28x10“ 12 
0.310 


Above  Front 


5.36x10“ 12 
0.022 


6.22x10“ 12 
0.011 


8.68x  10“12 
0.061 


2.36x  10-12 
0.118 


14.96 


1.24x  10“ 12 
0.022 


1.03x10“" 

0.039 


1.41x10“" 

0.668 


11.73 


1.12x10“" 

0.013 


6.33x  10“ 12 
0.121 


5.9  lx  10“ 12 
0.109 


12.79 


1.08x10“" 

0.015 


1.255x10  " 
0.082 


12.43 


Note:  the  numbers  below  normalized  dose  are  fractional  standard  deviation. 


Table  XI 


Normalized  Dose  Results  -  Series  #3 


Organ  # 


Muscle 


Brain 


ISOTROPIC 


6.94x  10“12 
0.046 


3.54x  10  12 
0.078 


4.02x10“ 12 
0.050 


6.90x  10“12 
0.172 


1.55x10“ 12 
0.401 


5.53 


Norm.  Dose  fGrav 


Rt.  Side 


1.07x10“" 

0.016 


5.51x10  12 
0.020 


6.43x  10“12 
0.016 


•  ni/cm-)}  from  given  source 


2.89x10“ 12 
0.096 


13.96 


Above 


3.46x  10“ 12 
0.017 


1.22x10“ 12 
0.028 


1.24x10“ 12 
0.020 


1.08x10 
0. 


2.90x10“ 13 
0.7-17 


Front 


1.45x10“" 

0.014 


7.47x  10“ 12 
0.023 


1.13x10“" 

0.015 


6.56x10  12 
0.096 


1  1 .92 


l.oOxlO 


8 


9.54x10“ 12 
0.018 


1.09x10“" 

0.013 


1.24x10“" 

.054 


12.00 


Note:  the  numbers  below  normalized  dose  are  fractional  standard  deviation 


Table  XII 


1 
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TO 

*4 


3 

I 


i 
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Normalized  Dose  Results  -  Series  jj=A 


Organ  # 

\orm.  Dose  [Grav/I  n  Inn  J)'|  from  given  source 

ISOTROPIC 

Rt.  Side 

Above 

Front 

Skin 

6.80x10“ 12 
0.036 

1.08x10"“ 

0.011 

3.50x10" 12 
0.013 

1. 47x10"“ 
0.011 

Bone 

3.35x10" 12 
0.053 

5.47xl0"12 

0.013 

1.24x10" 12 
0.018 

7.60xl0“12 

0.016 

Muscle 

4.14x  10“12 
0.037 

6.36xl0"12 

0.011 

1.24x10" 12 
0.015 

1.13x10"“ 

0.011 

Brain 

5.63x  10" 12 
0.133 

9.64xl0"12 

0.047 

l.llxlO"11 

0.024 

6.02x10" 12 
0.089 

GI  Tract 

1.66xl0"12 

0.244 

2.49xl0"12 

0.073 

1.16xl0~u 

0.612 

6.49x10“ 12 
0.073 

RSBIilH 

10.55 

26.89 

Back 


1. 22  x  10~ 11 
0.040 


6.67x  10~12 
0.064 


23.83 


Note:  the  numbers  below  normalized  dose  are  fractional  standard  deviation. 

In  all  four  cases,  the  normalized  doses  fluctuate  between  higher  and  lower 
values  with  respect  to  the  test  case  (Table  VI).  All  of  the  values  (except  gastro¬ 
intestinal  tract  from  above,  described  below)  are  within  acceptable  statistical  lim¬ 
its  of  the  test  case  values  showing  the  random  nature  of  a  Monte  Carlo  code. 
Although  the  standard  deviations  are  not  significantly  lower  compared  with  the 
test  case,  the  CPU  time  is  dramatically  reduced  in  series  2  and  3  (up  to  8oclj 
reduction).  The  very  slight  reduction  in  variance  resulting  from  the  additional 
neutrons  analyzed  in  the  fourth  series  does  not  justify  the  extra  CPU  time 
required  to  run  the  series.  The  solution  for  the  gastro-intestinal  tract  from  neu¬ 
tron  started  above  is  suspect  in  all  four  series  and  the  test  case.  Neutrons  started 
from  above  the  phantom  have  many  mean-free  paths  to  travel  before  reaching  the 
gastro-intestinal  tract.  Since  a  large  majority  of  the  neutrons  either  scatter  out  of 
the  phantom  or  are  statistically  removed  before  reaching  the  tract,  poor  statistics 
are  the  result  for  this  organ.  To  overcome  this,  many  more  neutrons  need  to  be 
started.  Unfortunately  the  computer  system  used  does  not  have  the  buffer  space 
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needed  to  make  such  an  accommodation. 

Sample  Calculation. 

A  sample  applications  problem  demonstrates  the  use  of  MORSE  phantom 
results.  Consider  a  10  kton  tactical  weapon  detonation  at  sea  level  on  solid,  flat 
ground.  The  problem  is:  "what  total  body  dose  does  a  soldier  1.5  km  from  ground 
zero  receive  if  the  blast  is  to  the  soldier’s  back?"  The  assumptions  required  are: 

1.  the  weapon  is  a  plutonium  fission  device, 

2.  neutrons  do  not  scatter  from  the  surface  of  the  earth, 

3.  the  atmosphere  is  homogeneous, 

4.  the  neutrons  are  mono-directional  when  they  reach  the  soldier, 

5.  the  neutrons  reaching  the  soldier  have  an  energy  distribution  as 
described  in  Table  V, 

6.  neutrons  have  an  average  quality  factor  per  energy  group  of  10  (12:373). 

A  47 rr2  normalized  fluence  per  source  neutron  of  0.02  (3:  Ch.9,18)  results 
from  a  1.5  km  target  range  at  sea  level.  The  number  of  neutrons  escaping  the 
device  is  1.1  X  1023  neutrons/kton  (3:Ch  9,2).  The  resulting  neutron  fluence 
incident  on  the  soldier  is  9.2  X  1010  n/cm2.  Figure  6  is  a  graphical  representation 
of  Table  VI,  the  original  analysis  results.  Reading  normalized  absorbed  dose 
values  in  the  -Y  direction  for  each  organ  off  the  graph  and  multiplying  by  the 
incident  fluence  results  in  an  absorbed  dose  of  approximately  5  Gray-tissue  (500 
Rad-tissue).  To  determine  an  equivalent  dose,  the  absorbed  dose  is  multiplied  by 
the  quality  factor  of  10.  The  equivalent  dose  is  50  Sievert  (5000  REM),  which 
would  render  the  average  human  incapacitated  within  approximately  30  minutes 


Figure  6.  Normalized  Absorbed  Dose  for  Test  Case 


IV.  Conclusions  and  Recommendations 

Conclusions 

The  goal  of  this  study  is  to  provide  a  fast,  convenient,  and  generic  method  of 
estimating  absorbed  neutron  radiation  dose  to  human  organs.  Modification  of  the 
ORNL  MORSE  code  proves  to  be  a  useful  technique  for  achieving  this  goal.  Con¬ 
clusions  about  the  user-ready  phantom  version  of  MORSE  are  as  follows: 

1)  The  most  efficient  MORSE  analysis,  to  date,  is  use  of  a  10  microsecond 
time  kill.  Time  kill  reduces  analysis  time  by  up  to  85%  without  degrad¬ 
ing  absorbed  dose  results  or  sample  variance. 

2)  Through  simple  changes  to  the  input  data  file,  normalized  absorbed  neu¬ 
tron  dose  estimates  are  easily  achievable  for  a  wide  variety  of  neutron 
environments.  Each  analysis  requires  under  15  CPU  minutes. 

3)  Although  this  version  of  MORSE  is  functional,  additional  modifications 
are  foreseeable  to  produce  a  more  efficient  and  meaningful  code.  The 
next  section  contains  recommendations  for  further  modifications. 

Recommendations 

The  following  six  recommendations  could  produce  a  much  more  efficient  and 
generic  phantom  code. 

Determine  a  method  of  increasing  the  number  of  allowable  neutrons  per 
batch.  This  could  improve  the  variance  in  the  smaller  organ  estimates. 

Follow  a  systematic  approach  to  reduce  the  sample  variance  using  split¬ 
ting,  Russian  roulette,  and  path  length  stretching.  Path  length  stretch¬ 
ing  is  not  attempted  in  this  project,  because  of  shortage  of  time,  but 
may  benefit  future  modifications. 

Implement  37  group  (or  less)  quality  factors  to  produce  effective  dose  in 
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Include  additional  phantom  geometries  such  as  kneeling  and  sitting  posi¬ 
tions.  The  additional  geometries  more  closely  resemble  common  battle 
situations  (i.e.  tank  driver  or  pilot). 

Include  more  organs  of  interest  if  somatic  radiation  effects  are  of  con¬ 


cern. 


Include  dose  from  internally  produced  gamma-rays  to  more  closely  esti¬ 
mate  actual  dose. 
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Appendix  A:  Example  of  MORSE  Input  Data  File 
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Appendix  A  is  a  listing  of  the  input  data  file  (excluding  cross  sections) 
required  by  MORSE.  The  format  of  the  data  file  is  given  in  the  MORSE  Applica¬ 
tions  Guide  (5).  The  file  is  set  up  to  produce  a  fission  neutron  energy  spectrum 
incident  on  a  three-dimensional  phantom. 
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MORSE,  PHANTOM  GEOMETRY,  37  n  groups 
1200  1200  20  1  37  0  37 


energy  spect. )  -  TEST  FILE 
0  0  80.  3  0 


p>l 

0 

37 

0  01 

.0  0 

.00001  0 

'.  0 

i.O  4380. 

IlA  * 

0.0 

0.0 

0.0 

0.0 

1.0 

1.0 

1.0 

.v.: 

0.0000 

+0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

Wjm  **  r 

0.0000 

+0 

3.8400  -3 

3.5000  -3 

5.3900  -3 

7.3500  -3 

1.8400  -2 

3.2500  -3 

■ 

8.4700 

-3 

5.5000  -2 

3.2400  -2 

1.0600  -2 

9.7200  -2 

1.4700  -1 

2.1600  -1 

1.5000 

-1 

1.9300  -2 

1.2100  -1 

5.7300  -2 

6.0000  -3 

2.4000  -2 

1.4400  -2 

rK 

0.0000 

+0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

0.0000  +0 

k* 

O.OCOO 

+0 

0.0000  +0 

1.9600 

+7 

1.6900  +7 

1.4900  +7 

1.4200  ->-'7 

1.3800  +7 

1.2800  +7 

1.2200  +7 

Ijfj 

1.1100 

+7 

1.0000  +7 

9.0500  +6 

8.1900  ^6 

7.4100  +6 

6.3800  +6 

4.9700  +6 

Ih 

4.7200 

+6 

4.0700  +6 

3.0100  +6 

2.3900  +6 

2.3100  +6 

1.8300  +6 

1.1100  +6 

B 

5.5500 

+  5 

1.5800  +5 

1.1100  +5 

5.2500  +4 

2.4800  +4 

2.1900  +4 

1.0300  +4 

3.3500 

+3 

1.2300  +3 

5.8300  +2 

1.0100  +2 

2.9000  +1 

1.0700  +1 

3.0600  +0 

'jjj 

1.1300 

+0 

4.1400  -1 

A 

000035FA731A 

\ 

0 

0 

0  0 

0  6 

1 

x 

0 

0 

0  0 

m 

0 

0 

Geometry  of  BODY  ( ski 

.n,  bones. 

AND  MUSCLE) 

B 

SPH 

1 

0.0 

2.5 

86.5 

9.0 

$ 

SPH 

2 

0.0 

2.5 

86.5 

8.8 

S 

RCC 

3 

0.0 

-9.0 

80.0 

0.0 

9.0 

-1.0 

& 

7.5 

RCC 

4 

0.0 

-8.8 

80.0 

0.0 

8.8 

-1.0 

7.3 

11 

RCC 

5 

0.0 

0.0 

69.0 

0.0 

1.0 

12. 

w 

5.2 

S3 

RCC 

6 

0.0 

0.0 

69. 

0.0 

1.0 

12. 

5.0 
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7 

0.0 

0.0 

0.0 

0.0 

0.0 

70.0 

Sc 

0.0 

10.5 

0.0 

15.5 

0.0 

0.0 

REC 

8 

0.0 

0.0 

0.2 

0.0 

0.0 

69.6 

■d  a 

0.0 

10.3 

0.0 

15.3 

0.0 

0.0 

TRC 

9 

-25.0 

0.0 

23.0 

12.0 

0.0 

45.0 

3.0 

4.0 

TRC 

10 

-25.0 

0.0 

23.0 

12.0 

0.0 

45.0 

Pi"** 

2.8 

3.8 

[  JS 

TRC 

11 

25.0 

0.0 

23.0 

-12.0 

0.0 

45.0 

3.0 

4.0 

■R 

TRC 

12 

25.0 

0.0 

23.0 

-12.0 

0.0 

45.0 

JE 

2.8 

3.8 

TRC 

13 

-25.0 

0.0 

-17.0 

0.0 

0.0 

40.0 

qS 

2.0 

3.0 

5H 

TRC 

14 

-25.0 

0.0 

-16.8 

0.0 

0.0 

39.8 

1.8 

2.8 

KS 

TRC 

15 

25.0 

0.0 

-17.0 

0.0 

0.0 

40.0 

™ 

2.0 

3.0 

It 

TRC 

16 

25.0 

0.0 

-16.8 

0.0 

0.0 

39.8 

p© 

1.8 

2.8 

B 

TRC 

17 

-10.0 

0.0 

-40.0 

2.5 

0.0 

42.0 

K 

6.0 

8.5 

K 

TRC 

18 

-10.0 

0.0 

o 

o 

1 

2.5 

0.0 

41.8 

D 

5.8 

8.3 

B 1 

TRC 

19 

10.0 

0.0 

-40.0 

-2.5 

0.0 

42.0 

9 

6.0 

8.5 

S 

TRC 

20 

10.0 

0.0 

-40.0 

-2.5 

0.0 

41.8 

AT 

5.8 

8.3 

2» 

TRC 

21 

-10.0 

0.0 

-80.0 

0.0 

0.0 

41.0 

V 
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K 
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REC  61 

0.0 

0.0 

9 . 0 

0.0 

9 . 0 

5.0 

0.0 

5.0 

0.0 

9.0 

0 . 0 

0 . 0 

SPH  62 

0.0 

0.0 

10.0 

101.0 

RPP  63 

-120.0 

120.0 

-120.0 

12C.0 

120.0 

1 2 

0 . 0 

END 

SKN  OR 

+  1 

-4 

-6 

-2CR 

+  3 

-2 

-6 

-40R  -5 

-2 

-4 

-8 

-60R 

+7 

-6 

-10 

_  1  v 

-18 

-20 

-26 

-80R 

+  9 

-3 

-10OR 

+  11 

-3 

—  12 

OR 

+  13 

-140R 

+  15 

-160R 

+17 

-8 

-26 

-180R  -19 

-3 

-26 

-20OR 

+  21 

-22i 

OR  +23 

-240R 

-25 

-3 

—  j. 

-19 

-25 

SON  OR 

-30 

-2SCR 

-29 

-23 

-3CCR 

-31 

-CSC 

?s  -3  2 

OR 

+  33 

-34 

-350R 

+  3  6 

-37 

-38 

-39 

-40 

—  4  2 

-42 

-43 

-44 

-45 

-46 

-47 

-48 

-49 

-50 

OR 

+  510R 

+  520R 

+  530R 

+  540R 

-550R  +560R 

+  570R 

+58CR  +59 

MUS  OR 

+  2 

-27 

-29 

-310R 

+4 

-290R 

-8 

-31 

-32 

-33 

-36 

-51 

-610R 

-33 

-35 

-3  4CR 

-3  6 

-33 

-370R 

+36 

+  39 

-3~CR 

-36 

-40 

-37CR 

-3  6 

-41 

-370R 

+36 

-42 

-37CR 

-36 

+43 

-37 

-51CR  -36 

-44 

-37 

-510R 

+  36 

+  45 

-37 

-510R 

-3  6 

-4  6 

-37 

-51CR 

+  36 

+  47 

-37 

-510R 

-36 

-43 

-37 

-510R 

+  36 

+  49 

-370R 

+  36 

+  50 

-370R 

+  28 

-60 

OR 

+  300R 

+  34 

-6  ICR 

+37 

-31 

-32 

-510R 

+  6 

-31 

-29 

-270R 

+  26 

-8 

-18 

-20OR 

-10 

-8 

-56 

OR 

+12 

-8 

-570R 

+  14 

-580R  +16 

-590R 

+  18 

-3 

-26 

-520R 

+  20 

-8 

-26 

-53GR 

-22 

-54 

-52 

OR 

+24 

-55 

-53 

BRN 

+28 

-60 

GIT 

+61 

AIR 

+62 

-1 

-3 

-5 

-7 

-9 

-11 

-13 

-15 

-17 

-19 

-21 

-23 

-25 

ENV  . 

+63 

-62 

END 
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3  4 

5 
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6 

1  2 

3  4 

5 

6 

0 

PHANTOM  PROBLEM  37 n, 

for  11  elements  ,  RUN  3 

37  37 
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58 

31 

4  3 

3 

3  4 

2 

0 

3 

0  0 

0  0 

0 

0 

0  -10 

0 

0  0 

SAMBO  ANALYSIS  INPUT 

DATA: 

phantom  prcble 

m 

5  1 

1  0 

0 

1 

1  1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

NORMALIZED  ORGAN  DOSES 

GRAY-tissue/( s-n/cm* 

*2)  [DET :  1 , 3  , 

4,5]  or 

GRAY 

-BONE/ ( s-n/cm*  *  2 ) 

[DET 

’:2] 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

USER  DEFINED 

INPUT:  : 

SOURCE 

TYPE  AND  X-SEC 

TION 

MEDIA 

12  11 

1  3 

3340.  6710.  49570. 

1555. 

2120. 

3.53E6 

•5' 

ft 
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Appendix  B:  Final  Version  of  Modified  MORSE  Subroutines 


Appendix  B  contains  the  MORSE  driver  program  and  six  subroutines 
modified  to  analyze  the  human  phantom  model.  The  subroutines  were  either 
modified  or  completely  re-written  to  provide  a  phantom  specific  code.  The  subrou¬ 
tines  are: 

GTMED 

-  determines  type  of  material  within  a  zone  (bone,  tissue,  air) 

BAXKR 

-  bookkeeping  subroutine 

TRKBDR 

-  determines  dose  contributions  from  particles  crossing  a  boundary 
TRKCOL 

-  determines  dose  contributions  from  particles  undergoing  a  collision 
INS  COR 

-  reads  additional  input  data  specific  to  this  study 
SOURCE 

-  produces  a  neutron  source  incident  on  the  phantom  either  isotropically 
from  a  sphere  surrounding  the  phantom,  or  mono-directionally  from  the 
right  side,  above,  front,  or  back  of  the  phantom. 

Detailed  descriptions  of  all  MORSE  subroutines  are  in  the  Applications  Guide 
(5).  Instructions  on  the  use  of  MORSE  are  in  the  thesis  notebook. 


Modifications  by: 

Capt  Michael  G.  Archuleta 
GNE-88M 

modifications  done  Oct/Nov  1987  for  thesis  project 
to  determine  organ  dose  as  a  function  of  n  flux. 

This  version  does  not  determine  uncollided  fluence. 
Instead,  a  tracklength  estimater  is  used  to  determine 
fluence  and  is  called  for  collisions  and  boundary  cross! 
Fluence  is  subsequently  converted  to  normalized  dose  in 
GRAY/(n/cm**2)  of  tissue  or  bone. 

This  version  requires  the  PHANTOM.DAT  input  data  file, 
PHANTSRC.FOR  source  generator,  and  MED3.BIN  binary  cross 
section  data  file. 


*  THIS  IS  THE  MAIN  ROUTINE  ********************* 

* 

*  THE  FOLLOWING  CARD  DETERMINES  THE  SIZE  ALLOWED  FOR  BLANK  COMMON  * 
COMMON  NC; 50000) 

*  (REGION  SIZE  NEEDED  IS  ABOUT  150K  -  4* (SIZE  OF  BLANK  COMMON  IN  WO 

*  NOTE  -  THE  ORDER  OF  COMMONS  IN  THIS  ROUTINE  IS  IMPORTANT  AND  MUST 

*  POND  TO  THE  ORDER  USED  IN  DUMP  ROUTINES  SUCH  AS  HELP,  XSCHL? ,  AN 

* 

*  LABELLED  COMMONS  FOR  WALK  ROUTINES  **************** 
COMMON  /APOLLO/  AGSTRT, DDF, DEADWT ( 26 ) , ITCUT  ,  ITIN 

COMMON  /FISBNK/  MFISTP 
COMMON  /NUTRON/  NAME 

★ 

*  LABELLED  COMMONS  FOR  CROSS-SECTION  ROUTINES  *********** 
COMMON  /LOGS LG/  ISCCOG 

COMMON  /MEANS/  NM 
COMMON  /MOMENT/  NMOM 
COMMON  /QAL/  Q 
COMMON  /RESULT/  POINT 

★ 

*  LABELLED  COMMONS  FOR  GEOMETRY  INTERFACE  ROUTINES  ******** 
COMMON  /GECMC/  XTWO 

COMMON  /NORMAL/  UNORM 

* 

*  LABELLED  COMMONS  FOR  USER  ROUTINES  *************** 
COMMON  /PDET/  ND 

COMMON  /USER/  AGST 

★ 

*  COMMON  /DUMMY/  WILL  NOT  BE  FOUND  ELSEWHERE  IN  THE  PROGRAM  *  *  *  * 
COMMON  /DUMMY/  DUM 

* 

*  LABELLED  COMMONS  FOR  USER  INPUT  DATA  ********* 

COMMON  /'INSDAT/ISRC, VOLI ( 10 )  ,IXSECM(10) 

CHARACTER* 20,  NAM1 
CHARACTER* 20,  NAM2 
TYPE  *, ’  ’ 

TYPE  *,'*********  MORSE  Code,  HUMAN  PHANTOM  Problem  *********' 

TYPE  *,'  - - - — ->  WARNING  !!!  <- — ■ — ' 

TYPE  *, 'ABORT  if  mixed  x-secs  are  not  assigned  to  FOR010' 

TYPE  *, ‘  ' 

TYPE  *, 


ENTER  NAME  OF  INPUT  FILF' 


ACCEPT  100,  NAM1 
FORMAT (A20) 


nnnn 


TYPE  -*  ,  '  ENTER  NAME  OF  OUTPUT  FILE' 

ACCEPT  200,  NAM2 
200  FORMAT  (A20) 

OPEN( UNIT-1, NAME-NAM1, TYPE- 'OLD' ) 

OPEN ( UNIT- 2 , NAME-NAM2 , TYPE- ' NEW ' ) 

ITOUT  -  2 
ITIN  -  1 
NLFT-34999 
CALL  MORSE (NLFT) 

TYPE  jOU,  NAM 2 

300  FORMAT ( X, 'OUTPUT  FILE  IS  ' ,A20) 

STOP 

END 

£*,**■***•****•*******•»»•*  subroutine  GTMED*’'’"***”'”*********'’ 

SUBROUTINE  GTMED(MDGECM,MDXSEC) 

£***•*★*'■**■* 

C  **  This  version  of  GTMED  assigns  xsecticn  media  to  geometry  media 
C  **  as  specified  in  the  additional  user  input,  lines  in  the  input 
C  **  data  file 

COMMON  /INSDAT/ISRC, V0LI ( 10 ) , IXSECM( 10 ) 

IF(MDGECM.EQ. 0 . OR.MDGEOM. EQ. 1000 )GO  TO  90 
MDXSEC  -  IXSECM(MDGECM) 

RETURN 

90  MDXSEC  -  MDGECM 
RETURN 
END 

C  ************** FUNCTION  DIREC*  * *  +  ■*•*■**■*****■*■**■*** 

FUNCTION  DIREC(F) 

******** 

**  This  version  of  DIREC  is  a  dummy  routine  that  does 
**  absolutely  nothing. 

******** 

direc-1 , 

RETURN 
END 

-■.***■»*.»*■****■»» SUBROUTINE  BANKR*  * **»»**■*■***•*■****■* 

SUBROUTINE  B ANKR ( NBNKID ) 
c********* 

C  **  This  version  call  TRKBDR  from  BANKR ( 7 )  and  TRKCOL  from  BANKR ( 5 ) . 

£  i r******** 

COMMON  /APOLLO/  AGSTRT , DDF , DEADWT (  5 )  , ETA , ETATH , ETAUSD , UINP , VINP , 

1  WINP , WTSTRT, XSTRT, YSTRT, ZSTRT, TOUT, XTRA( 10 ) , 

2  10,  II, MEDIA, IADJM, IS3IAS , ISOUR, ITERS, ITIME, ITSTR, LOCWTS , LOCFWL, 

3  LOCEPR , LOCNSC , LOCFSN , MAXGP , MAXTIM , MEDALB , MGPREG , MXREG , NALB , 

4  NDEAD(  5 ) ,  NEWNM , NGEOM , NGPQT1 , NGPQT2 , NGPQT3 , NGPQTG , NGPQTN, NITS , 

5  NKCALC ,  NKILL ,  NLAST,  NMEM,  NMGP ,  NMOST ,  NMTG ,  NOLEAK ,  NORMF ,  NPAST , 

6  NPSCL( 13 ) , NQUIT, NSIGL , NSOUR , NSPLT , NSTRT , NXTRA ( 10 ) 

COMMON  /NUTRON/  NAME , NAMEX , IG , IGO , NMED , MEDOLD , NREG , U , V , W , UOLD , VOLD 
1  ,  WOLD ,  X, Y, Z , XOLD , YOLD , ZOLD , WATE , OLDWT , WTBC ,  BLZNT , BLZON , AGE , OLDAGE 
NBNK  -  NBNKID 
IF  (NBNK)  100,100,140 

100  NBNK  -  NBNK  +  5 

GO  TO  (104,103, 102,101) , NBNK 

101  CALL  STRUN 

C  CALL  HELP (4HSTRU, 1,1, 1,1) 

RETURN 

102  NBAT  -  NITS  -  ITERS 
NSAVE  -  NMEM 

CALL  STBTCH(NBAT) 


nnnnnnnn  noonn 


C  N3AT  IS  THE  3ATCH  NO.  LESS  ONE 
RETURN 

103  CALL  NBATCH ( NSAVE ) 

C  NSAVE  IS  THE  NO.  OF  PARTICLES  STARTED  IN  THE  LAST  BATCH 
RETURN 

104  CALL  NRUN ( NITS , NQUIT ) 

NITS  IS  THE  NO.  OF  BATCHES  COMPLETED  IN  THE  RUN  JUST  COMPLETED 
NQUIT  .GT.  1  IF  MORE  RUNS  REMAIN 

.EQ.  1  IF  THE  LAST  SCHEDULED  RUN  HAS  BEEN  COMPLETED 
IS  THE  NEGATIVE  OF  THE  NO.  OF  COMPLETE  RUNS,  WHEN  AN 
EXECUTION  TIME  KILL  OCCURS 


RETURN 

140  GO  TO  (1,2,3,4,5,6,7,8,9,10,11,12,13), NBNK 


NBNKID 

1 

3 


13 


COLL  TYPE 
SOURCE 
FISSION 
REAL  CCLI 
BDRYX 
K-— Cur 
R  R  KILL 
GAMLCST 


1 

2 

3 

4 

5 


3ANKR  CALL 
YES  ( MSOUR ) 
YES  ( FPROB ) 
YES  (MORSE) 
YES  (NXTCOL) 
NO  (MORSE) 
(TESTW) 
(GSTCRE) 


NBNKID 

2 

i 


COLL  TYPE 
SPLIT 
G  AMGEN 


3ANKR  C 
NO  (TE 
YES  (GS 


STW) 

TORE 


vrc 


NO 

NO 


R  R  SURV 


RETURN 
RETURN 
RETURN 
RETURN 
CALL  TRKCOL 
RETURN 

6  RETURN 

7  CALL  TRKBDR 
RETURN 

8  RETURN 

9  RETURN 

10  RETURN 

11  RETURN 

12  RETURN 

13  RETURN 
ENT) 

£♦*■*•***'*■*■*•*•*■**•*'**•*'*•*'*•*•**■« 
SUBROUTINE  TRKBDR 


SUBROUTINE  TRKBDR  ******* 


4r  x  x  x  ■w  1 


Qx  *  *  *  *  x  x  x  x  *  X  *  *  *  x  *  *  *  x  *  *  x  *  *  *  *  x  *  *  *  *  *  *  *  *  *  *  x  x  x  x  x  *  x  *  *  x  i 


tXXXXXXXXXXXf 


this  version  determines  flux  as  tracklength  divided  by  detector 
volume  for  all  boundary  crossing  out  of  body  organs, 
the  variables  used  are: 


AREA 

CON 

ISRC 

KERMAB 

KERMAT 


-  surface  area  of  neutron  source  region  [cm** 2] 

-  dose  contribution  of  n  in  MEDOLD  [GRAY/(n/cm**2) ] 

-  type  of  source  defined  in  input  data  file 

-  bone  KERMA  factor  [RAD-bone/(n/cm**2) ] 

”  tissue  KERMA  factor  [RAD-tiss/(n/cm**2) ] 

MEDOLD  -  media  (organ)  n  is  exitting 
MXREG  -  highest  region  number  (corresponds  to  air) 

TRK  -  tracklength  from  previous  collision/boundary 
crossing  to  current  boundary  crossing  [cm] 

VOL  -  volume  of  media  (organ)  n  is  exitting 

Q*  *****  X  ****************  x  ***************  X  *  *  X  X  *******  X  *****  X  *  X  ******  X 

COMMON  /APOLLO/  AGSTRT , DDF , DEADWT ( 5 ) , ETA , ETATH , ETAUSD , UINP , VINP , 

1  WINP , WTSTRT , XSTRT , YSTRT , ZSTRT , TCUT , XTRA  ( 10 )  , 

2  10,  II, MEDIA, IADJM, ISBIAS , ISOUR , ITERS , ITIME, ITSTR , LOCWTS , LCCFWL , 

3  LCCEPR , LCCNSC , LCCFSN , MAXG? , MAXTIM , MEDAL3 , MGPREG , MXREG , NAL3 , 

4  NDEAD (  5  )  ,  NEWNM,  NGECM,  NGPQT1 ,  NGPQT2  ,  NGPQT3  ,  NGPQTG  ,  NGPQTN ,  NITS  , 

5  NKCALC , NKILL, NLAST, NMEM , NMG? , NMCST ,  NMTG , NOLEAK , NORMF , NPAST , 


nnnnnnnn 


6  N?SCL(  13  )  ,  NQUIT,  NSIGL,  NSCUR ,  NSPLT ,  NSTRT  ,  NXTRA(  10  ; 

COMMON  /NUTRON/  NAME , NAMEX , IG . IGC , NMED , MEDOLD , NREG , U , V , W , UCLD - YCLD 
1  ,  WOLD ,  X,  Y, Z ,  XOLD , YOLD , ZOLD , WATE , OLDWT , WT3C , 3LZNT , 3LZCN ,  AGE , CLDAGE 
COMMON  /INSDAT/ISRC , VOLI (  10  )  , IXSECM( 10) 

REAL  KERMAT (37), KERMAB (37), AREA ( 5 ) 

DATA  AREA/62800. ,4580. ,1560. ,10600. ,10600./ 

DATA  KERMAT/. 7 225E-08, .6858E-08, .6642E-08, .S528E-08, .6380E-08, 


* . 6188E-08, 

6118E-08 , . 5357E-08 , 

. 2  o  94E-0 8 , 

. 3393E-08, 

. 5284E-C8, 

* . 5039E-08 , 

4723E-08 , . 4479E-08 , 

. 4412E-09 , 

.4193E-08, .3560E-08, 

*  .  3264E-08, 

.  3200E-08, .276IE-08, 

. 2097E-08 , 

.1302E-08, 

. 8010E-29 , 

* . 5535E-09, 

.  3148E-09, . 2145E-Q9, 

. 14~2E-09 , 

. 6269E-1C , . 2211E-1: . 

* . 91S7E-11, 

* . 3689E-11, 

.3I04E-11, . 10  27E-1I , 
1515E-10/ 

.  1 3  •  E-  -  w  / 

. 1341E-I1, . 2251E-11 , 

DATA  KERMAB/. 5511E-08, . 5196E-C 

8 , . 5016E-C 

8, .4929E-08, . 4824F-: 

* . 4680E-08 , 

4608E-08, . 4396E-08 , 

. 4272E-08 , 

. 3991E-08, 

. 3911E-33 , 

* . 3663E-08  , 

3405E-08, . 3199E-08 . 

. 3162E-03  , 

. 2995E-08  ,  . 24893-33 , 

* . 2256E-Q8 , 

2214E-08, . 1906E-C3 , 

'  ^  p 

. 5  34  2E-39.  . 5 4  ~  3  E -  2  ? . 

* . 3308E-09 , 

2144E-09, . 14  6  IE— 39  , 

. 12C2E-09 , 

. 42“6E-1C , . 1321E-13 , 

* . 6449E-11, 

* . 6100E-11, 

2407E-11, . 1196E-11, 
2474E-10/ 

.  1415S-11 

. 2223E  11, 

.  3  7*4  6  L  - 1 1 , 

C  *  ♦  *  check  for  neutron  coming  from  VOID,  EN7,  or  AIR  *  T  * 

IF ( MEDOLD . EQ . MXREG . OR . MEDOLD . EQ . 0 ) RETURN 
C  *****  add  crack  length  to  zone  detector  *  *  *’,*****•*»•****’'**’' 

C  FLUENCE  [CON]  -TRK/VOL 

C  Norm.  DOSE  [rad/(n/cm**2)  ]  -  flue.nce  *  kerma  factor*area  * 

C  Gray/ ( n/cm* * 2 )  -  fluence  *  kerma  factor  *area/100.  * 

TRK  -  WATE  *  SQRT( (X-XOLD) **2  -  (Y-YCLD)**2  *  (Z-ZCLD)**2) 

VOL  -  VOLI (MEDOLD) 

IF ( MEDOLD. EQ. 2) CON  -  TRK  *  KERMA3 ( IG )  *  AREA( ISRC )/(VOL*100 . ) 

IF ( MEDOLD. NE. 2) CON  -  TRK  *  KERMAT ( IG )  *  AREA( ISRC )/( VCL *100 . ) 

C  if  source  is  isotropic,  correct  for  2pi  Drcaucticn 
IF( ISRC. EQ.l) CON  -  CON/2. 

CALL  FLUXST ( MEDOLD, IG, CON, 0.0, 0.0,0) 

C  *  *  SWITCH  -  0  —  STORE  IN  ALL  RELEVANT  ARRAYS  EXCEPT  UD 

RETURN 
END 

£******»*■*■»*, ■****•***•.*■*  SUBROUTINE  TRKCOL  **»*''***»**■***»«*** 

SUBROUTINE  TRKCOL 

-**-**-**-*-****'*ir-**-*ir*-**'*-*ir*-***i'ir-**T-*x-*x  +  xxir**  +  xxx*-*-Kxx*f*-xicKf*iri'x*xxic-xx-*x 

this  version  determines  flux,  at  each  collision  site, 
from  tracklength  divided  by  detector 

volume  and  in  used  with  TRKBDR  (called  from  BANKR ( 7 ) ) 

Variables  used  are  identical  to  TRKBDR  variables  except: 

NMED  -  current  media  (organ)  n  is  in. 

IGO  -  n  energy  before  collision 

★  *****★**★★**★**★****★★★★*★★**■*★*★**•*•*■**■*★***■**★*■*★*■*★**•**★★*★***■****'* 
COMMON  /APOLLO/  AGSTRT, DDF, DEADWT( 5 ) , ETA, ETATH , ETAUSD , UINP , VINP , 

1  WINP , WTSTRT, XSTRT, YSTRT, ZSTRT, TOUT, XTRA( 10 )  , 

2  10, II, MEDIA, IADJM, IS3IAS , ISOUR , ITERS , I TIME , ITSTR , LOCWTS , LCCFWL , 

3  LOCEPR , LOCNSC , LOCFSN , MAXGP , MAXTIM , MED ALB , MGPREG , MXREG , NAL3 , 

4  NDEAD( 5 ) , NEWNM, NGEOM, NGPQT1 , NGPQT2 , NGPQT3 , NGPQTG , NGPQTN, NITS , 

5  NKCALC , NK ILL , NLAST , NMEM , NMG  P , NMCST , NMTG , NOLEAK , NORMF , N  P AS  T , 

6  NPSCL( 13 ) , NQUIT , NSIGL , NSOUR , NSPLT , NSTRT , NXTRA(  10 ) 

COMMON  /NUTRON/  NAME , NAMEX , IG , IGO , NMED , MEDOLD ,  NREG , U, V, W, UOLD , VOLD 
1  ,  WOLD ,  X, Y, Z , XOLD , YOLD , ZOLD , WATE , OLDWT , WTBC ,  SLZNT , BLZCN , AGE , OLDAGE 
COMMON  /INSDAT/ISRC, VOLI (10) , IXSECM(IO) 

REAL  KERMAT (37), KERMAB (37), AREA  ( 5 ) 

DATA  AREA/6  2800. , 4  580 . ,  1560 . ,  106C0 . , 1C6C0  ./ 

DATA  KERMAT/. 7225E-08, . 6858E-Q8, . 6642E-08, . 6528E-08, .6380E-Q8, 


u  -  M 


♦.6I88E-08,  . 6 113  E- 08 ,  .5857E- 


* . 3 264 E- 08 


*  .  3689E-11 , . 1515E-10/ 


08  ,  . 

3694E-CS  ,  . 

5  3  9  3  E- 0 3 , 

.  5294E-08 , 

08  ,  . 

4412S-38  ,  . 

4193E-08, 

. 356CE-08, 

08  ,  . 

2097S-08  ,  . 

1302E-08, 

. 8010E-09 , 

09,  . 

14720-09  ,  . 

6269E-10, 

. 2211E-10, 

11,  . 

9187S-12,  . 

134 IE- 1 1 , 

. 2251E-11, 

E-08 

, . 3016E-08 

, . 4929E-08 , . 4824E-08 

08,  . 

4272E-08, . 

3991E-08, 

. 3911E-C3, 

08,  . 

3162S-08,  . 

2995E-08, 

. 2489E-03 , 

08,  . 

1445E-08,  . 

3942E-09, 

. 547CE-39, 

09,  . 

1C02E-09 , 

4276E-10 . 

. 1521E-10, 

1  1 

"  4 1 5  E- ' '  , . 

2223E-11, 

. 3746 E- 11, 

AREA(  ISRC)// 


*  .  6100E-11 , .2474E-10/ 

C  *  ’  *  check  for  ENV  OR  AIR 

IF ( NMED . EQ . MXREG . OR . NMED . EQ . 0 ) RETURN 
C  *  *  *  calculate  fluence  estimate 

TRK  -  WTSC  *  SORT ( ( X-XCLD )  *  *  2  -  (Y-YCLD;**2  -  ( Z-ZGLD) * '2 ) 

VOL  -  VOLI( NMED) 

IF ( NMED . EQ . 2 ) CON  -  TRX  *  XERMA3 ( IGO )  *  AREA( ISRC )/( VCL* ICO . ) 

IF ( NMED .  NE .  2 ) CON  -  TRX  *  XERMAT ( IGO )  *  AREA/ I5EC)/(YCLT120 .  } 

IF ( ISRC . EQ . 1 ) CON  -  CON/2.0 

CALL  FLUXST(NMED, IGO, CON, 0.0, 0.0,0) 

C  *  *  SWITCH  -0  —  STORE  IN  ALL  RELEVANT  ARRAYS  EXCEPT  UD 

RETURN 
END 

****************  *  *  SUB  ROUTINE  INSCOR*  **■***■*■*******'* 

SUBROUTINE  INSCOR 

Qt******imr****t?r**t*tt*****ir***»*'<ftTt»tTTT*tTt*t*tttt»»*t 

C  This  version  reads  in  4  lines  of  user  input  data  frcm  the  input 
C  data  file. 

C  LINE  1  -  dummy  comment  line 

C  LINE  2  -  The  type  of  source  incident  on  the  phantcm: 

C  1  -  isotropic 

C  2  -  source  plane  in  +X  direction 

C  3  -  source  plane  in  -X  direction 

C  4  -  source  plane  in  +Y  direction 

C  5  -  source  plane  in  -Y  direction 

C  LINE  3  -  assigns  cross-section  media  to  geometry  media 

C  LINE  4  -  defines  the  volume  of  each  region  (organ) 

Q*-**'***'*********r**-*Te***-***-****iirTr**'*irir-*ir-ir-«r*it'r**ir*'**'**-*-*Tr*** 

COMMON  /APOLLO/  AGSTRT , DDF , DEAEWT ( 5 ) , ETA,  ETATH , ETAUSD , UINP , VINP , 

1  WINP , WTSTRT , XSTRT , YSTRT , ZSTRT , TOUT , XTRA ( 10 ) , 

2  10, II, MEDIA, IAEVTM, ISBIAS , ISCUR , ITERS , ITIME , ITSTR , LOCWTS , LOCFWL , 

3  LOCEPR , LCCNSC , LOCFSN , MAXGP , MAXTIM , MEDALS , MG P REG , MXREG , NALB , 

4  NDEAD ( 5 ) , NEWNM , NGEOM, NGPQT1 , NGPQT2 , NGPQT3  ,  NGPQTG , NGPQTN , NITS  , 

5  NXCALC , NKILL , NLAST, NMEM , NMGP , NMOST , NMTG , NOLEAX , NORMF , NPAST , 

6  NPSCL ( 13 ) , NQUIT, NS IGL , NSOUR , NS ?LT , NSTRT ,  NXTRA  (10) 

COMMON  /INSDAT/ISRC ,  VOLI  ( 10  ) ,  I-JCSECM(IO) 

CHARACTER* 80  NSTRING 

CHARACTER* 15  TYPSRC 
READ( II, 10) NSTRING 

10  FORMAT (A80) 

READ(  II, 11) ISRC 

11  FORMAT  (12) 

IF( ISRC . EQ . 1 )TYR£RC  -  ’ISOTROPIC.' 

IF ( ISRC. EQ. 2) TYPSRC  -  'MONO-DIR  IN  ^X.  ' 

IF ( ISRC . EQ . 3 ) TYPSRC  -  'MONO-DIR  IN  -Z.' 

IF( ISRC. EQ. 4) TYPSRC  -  'MONO-DIR  IN  -Y. ' 

IF( ISRC. EQ. 5) TYPSRC  -  'MONO-DIR  IN  -Y.  ' 

WRITE ( 10, 21)TYPSRC 


Sup-10 


S'. 

*V> 


21  FORMAT (/5X,  ’  """""ADDITIONAL  DATA  READ  THROUGH  INSCOR*****’ 
*8X,'The  source  being  used  is  ’,A15) 

READ (11,12)  (IXSECM(I)  , I-1,MXREG) 

12  FORMAT ( 1513) 

WRITE<I0,22) (IXSECM(I) ,I-1,MXREG) 

22  FORMAT (/8X, 'THE  CROSS-SECTION  MEDIA  OF  EACH  REGION  ARE: '/ 

*  12X, 1513 ) 

READ( II , 13 ) ( VOEI ( I ) , 1-1 ,MXREG ) 

13  FORMAT (10E8. 3) 

WRITE ( 10,23) (VOLI ( I )  , I”1 ,MXREG) 

23  FORMAT (/8X, 'THE  VOLUMES  OF  EACH  REGION  ARE:'/ 

*  12X, 5E8 . 3/12X, 5E8 . 3 ) 

RETURN 

END 


'//, 


<• 


non  oonnonnnn 


SUBROUTINE  SOURCE ( IG,U,  V,  W,  X ,  Y,  Z ,  WATS , MED , AG , I SOUR, ITSTR, NGPCT3 - 
1  DDF, ISBIAS, NMTG) 

COMMON  /USER/  DUM( 9 ) , 10 , II, IDUM( 12 ) 

COMMON  WTS(l) 

COMMON  /INSDAT/ISRC , VOLI ( 10 ) , XSEC ( 10 ) 

This  version  produces  a  source  defined  in  the  user  input  section 
of  the  input  data  file.  The  possible  source  types  are: 

ISRC  -  1  -  isotropic 

2  -  source  plane  in  +X  direction 

3  -  source  plane  in  -Z  direction 

4  -  source  plane  in  -Y  direction 

5  -  source  plane  in  -Y  direction 


DATA  I CALL/1/ 

IF  (ICALL)  10,10,5 
I CALL  -  0 
WRITE  (10,1000) 

FORMAT  ( '  YOU  ARE  USING  THE  DEFAULT  VERSION  OF  SOURCE  WHICH  SETS  W 
1ATE  TO  DDF  AND  PROVIDES  AN  ENERGY  IG.') 

10  IF( ISOUR) 15,15, 60 
WATE-DDF 

IF  (ISBIAS)  20,20,25 
NWT  -  2* NMTG 
GO  TO  30 
25  NWT  -  3*NMTG 

30  R  -  FLTRNF ( 0 ) 

DO  35  1-1 , NGPQT3 
IF  (R  -  WTS( I+NWT) )  40,40,35 
35  CONTINUE 
40  IG-I 

IF  (ISBIAS)  60,60,45 
45  IF  (1-1)  60,50,55 

50  WATE  -  WATE*WTS ( 2*NMTG+1 )/WTS ( 3  *NMTG+1 ) 

GO  TO  60 

55  WATE  -  WATE*(WTS(2*NMTG+I)-WTS(2*NMTG+I-1) )/(WTS( 3*NMTG-I)-WTS( 3*N 
1MTG+I-1) ) 

*■*★*★***■*★■*★*★*■*★*•*■*★*★★*■*★*•*•*■* 

User  defined  source  description 

★  ★★★x-**-**********-************** 


60  GO  TO  (61,62,63,64,65) , ISRC 

C  **********  isotropic  source  *********** 

61  CALL  GTISO(XX,YY,ZZ) 

X  -  100.  *  XX 

Y  -  100.  *  YY 

Z  -  100.  *  ZZ  +  10. 

71  CALL  GTISO(U, V, W) 

DOTP  -  X*U  +  Y*V  +  Z*W 
IF ( DOTP . GE .0.0) GO  TO  71 
RETURN 

C  *********  source  plane  (+X)  *********** 

62  Y  -  26.  *  FLTRNF ( 0 )  -  13. 

Z  -  176.*  FLTRNF ( 0 )  -  80. 

X  - 30. 

U  -  1.0 

V  -  0.0 

w  -o.o 

RETURN 

C  *********  source  plane  (~Z)  *********** 

63  X  -  60.  *  FLTRNF ( 0 )  -  30. 


L3  . 


Y 

26.  * 

FLTRNF ( 0 ) 

2 

- 

96. 

U 

- 

0.0 

V 

- 

0.0 

w 

-1.0 

RETURN 

c 

★  ★**★**•** 

source 

plane  (+Y) 

64  X 

- 

60.  * 

FLTRNF ( 0 ) 

Z 

- 

176.  * 

FLTRNF ( 0 ) 

Y 

- 

-13. 

U 

- 

0.0 

V 

- 

1.0 

w 

- 

0.0 

RETURN 

c 

source 

plane  (-Y) 

65  X 

- 

60.  * 

FLTRNF ( 0 ) 

Z 

176.* 

FLTRNF ( 0 ) 

y.‘ 
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are  limited  to  organs  capable  of  incapacitiating  a  human.  The  code 
allows  for  five  different  source  direction  configurations  that 
simulate  neutrons,  of  any  specified  energy  distribution,  incident  on 
the  phantan. 

MORSE  analysis  of  a  fission  neutron  spectrum  on  the  phantcm  pro¬ 
duces  absorbed  dose  estimates  comparable  with  Japanese  atomic  bcmb 
survivor  dose  estimates  by  Scientific  Applications  International 
Corporation.  The  analysis  of  24,000  source  neutrons  requires  less 
than  15  central  processing  unit  minutes  on  a  VAX  11/780  computer 
(VMS  operating  system).  Although  the  code  is  currently  usable,  add¬ 
itional  phantcm  model  orientations,  ene  rgy -dependent  quality  fac¬ 
tors,  and  implementation  of  secondary  gamma-ray  dose  estimation 
could  greatly  improve  the  flexibility  and  usefulness. 


